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□ Filtration of VHF Sources 

□ Flash-Typing 

□ Transformations, Spatial 
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□ Channel Segment Creation 
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□ NOx Computation 

> Lab [Wanget al., 1998] 

> Theory [Coo ray et al. 2009] 
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Initial Processing 

1. McCaul VHF Source Clustering Algorithm 
S Ingests LMA data 

S Clusters (chi-square < 2) vhf sources into "Things" 

o "Things" are: Flashes, Small (Non-lightning) Discharges, Noise 

2. LNOM "Thing" Filtration 

S Range Filter (remove "thing" if completely outside cylinder) 

Number Filter (remove "thing" if <20 vhf sources... e.g., singletons) 

S What remains: legitimate flashes that are at least partially in cylinder 

3. LNOM VHF Source Filtration 

S Power Filter (remove vhf source if power < 1 dbW) 

S Altitude Filter (remove vhf source if at sfc or >=20km) 

4. LNOM Flash-Type Categorization 

S Deemed a CG if a NLDN CG detection is within (100 ms, 10km), and 
lowest altitude vhf is below N-region height. 

S But, deemed ambiguous if peak current is between 10 - 20 kA. 

5. LNOM Processing of VHF Sources in Legitimate Flashes 
S Coordinate Transformation 

S Spatial Averaging 
S Sorting 


Channel Length Algorithm 


o Begin with n averaging points (n = 5 in example below) 
o 1 st Iteration: 

■ Start at the highest altitude point. Define it to be "on the 
channel" and all other points to be "free" (off channel). 

■ Draw line from starting point to closest free point. This 
is the first channel "section". 

■ Now there are 2 points on the channel, and n-2 free points 


Starting point 


Channel Length Algorithm (cont.) 


o 2 nd Iteration: 

■ Find closest free point to 1 st channel point 

■ Find closest free point to 2 nd channel point 

■ Pick min of the mins (i.e., draw line from a free point to a 
channel point that is the minimum distance). 

o Continue with more iterations until no more free points 



Starting point 


Geodetic Height (km) Geodetic Height (km) 


Sample Cloud Flashes 

cloud L= 16851. 886 #aourcea=18 cloud L=26500.066 #aourcea=18 



cloud L= 1741 0.784 j|!aourcea=17 




cloud L= 14380.573 j|laourcea=12 
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Figure 1. Sample channel length, L, computed for 4 cloud flashes. Number of sources shown in upper 
right is the number of spatially averaged points (black dots). Green dots are the VHF sources that meet 
the 2 dBW power threshold, and the blue dots are those that do not. 
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Sample Ground Flashes 




L= 66087.061 #sourcea=20 




L=44866.474 #aources=39 




ground 


L=43004.363 


#aources=29 


ground 


L=41 106.261 #*ource»-3i 




Figure 2. Sample channel length computations for 4 ground flashes. 
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Channel length (August 2006) 
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STze (stots) = 8053 
STze (sho^n) = 7155 
Average = 51,438659 
Std Dev = 54,374795 
Max = 1350,7276 
MedFan = 36,204969 
Min = 1,8396797 



Computing egment Altitude Distribution 
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>Chop up Channel into 10-m Segments 

>Tally segments in each 100 m layer 
(only within cylinder) 
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x profile (August 2006) 
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Figure 9. Maximum increase in surface O s due to lightning NOx (ppb). 
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Figure 10 Maximum increase in boundary layer O s clue ro HghtningNOx tppb). 
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Figure 11. Maximum domain-wide increase in O? due to lightningNOx (ppb). 



NOx Profiles for all Months in 2007 (Latest LNOM Version) 
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Desired Future Application 



Ground Flash Fraction 
Retrieval Algorithm 


GLM = Geostationary Lightning Mapper (on GOES-R) 
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Air Quality Models (e.g. cmaq> 
Global Chemistry/Climate Models 


(e.g. GISS Model E, Geos Chem) 









Additional Current/Future Work 


We have now completed all months in 2005, 2006, and 2007. 
Working on 2008, 2009, and 2010. 

- My work over the past several months has been to complete flash 
clustering analysis for 2005-2010. 

- Currently working on DCLMA analysis for 2007. 

LNOM results being made available to multiple research partners: 

- NASA GSFC 

- UAH 

- Georgia Tech 

- Harvard 

- NASAGISS 

LNOM results will be placed on Global Hydrology Resource Center 
(GHRC) website to better serve customers. 


